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Abstract 

The critical behaviour of a system of two coupled scalar fields in three dimensions 
is studied within the formalism of the effective average action. The fixed points of 
the system are identified and the crossover between them is described in detail. 
Besides the universal critical behaviour, the flow of the coarse grained free energy 
also describes the approach to scaling. 
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1. Introduction 



The properties of multicomponent statistical systems near the critical temperature of a 
second-order phase transition are often characterized by a crossover behaviour. Crossover 
typically occurs if two or more fixed points are relevant for the description of scaling. Ex- 
amples can be found in physical systems with small impurities or a small coupling between 
two degrees of freedom. A well developed formalism for the description of the universal 
features of critical crossover is based on the renormalization group equations [0-0. Nev- 
ertheless, the practical situation in a given experiment may be quite complicated. The 
correlation length is large but not infinite in real life and, as a consequence, the scaling 
properties are not perfect. The temperature dependence of order parameters, susceptibil- 
ities, etc., often reflects an interplay between the approach to scaling and crossover. The 
first is determined by the flow of couplings towards a fixed point, whereas the second cor- 
responds to the running of couplings from one fixed point to another. Standard methods 
employed for the study of critical phenomena, such as for example the e-expansion, are 
often insufficient for the quantitative description of such complex situations, which are in 
addition not completely governed by universal scaling behaviour. Numerical simulations 
are also difficult in this context, as typically several different length scales are involved. 

We propose here the use of a field theoretical description in terms of the effective 
average action IV This corresponds to a coarse grained free energy where fluctuations 
with momenta larger than k are integrated out. The flow of I\ with a change of the 
infrared cutoff scale k is described by an exact non-perturbative flow equation. This can 
be solved approximately by appropriate truncations of the most general form of IV At 
a given scale k only modes with momenta around k contribute to the flow. We specify 
the microscopic physics by the "initial value" of IV at some high momentum scale A. For 
a given experimental situation this step involves the rewriting of the problem in terms 
of continuous variables. This part of the problem is not universal and involves only the 
short distance modes. It is not treated in the present paper. We only mention here that, 
for lattice models, it is advantageous to take A somewhat smaller than the inverse lattice 
spacing, so that deviations of the form of T\ from translational and rotational invariance 
are not too large. Our formalism can be extended to situations without translational and 
rotational symmetry. It is, however, much simpler for translationally and rotationally 
invariant situations where standard field theory can be used. In the vicinity of the critical 
temperature T cr of a second-order phase transition one of the couplings parametrizing IV 
- in our case 5k,\ - can be used to describe the distance from the critical hypersurface in 
the space of couplings, i.e. Sk\ ~ T — T cr . The order parameter, correlation lengths etc. 
can then be studied as functions of 5k^- For crossover situations there exists a second 
important coupling - in our case x\ - which determines to which fixed point the flow of 
the couplings is dominantly attracted. The parameter xa typically measures the amount 
of impurities or the size of a coupling between different components of a system. Critical 
quantities, such as the exponent v which governs the dependence of the correlation length 
on {y = — <91n£/<91n |5ka|), now depend on the ratio x\/\6k\\. 

We determine the long distance physics by solving numerically the flow of r&, starting 
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at short distances with Ta and moving towards k = 0. For k = the effective average 
action equals the free energy expressed as a functional of the fields. In the language of 
field theory, To is the generating functional of the 1PI Green functions, and all physically 
relevant information can be extracted from it. The method is formulated in arbitrary 
dimensions and can, in particular, be applied directly to three-dimensional statistical 
systems. Consequently, there is no restriction to the scaling situation for innnitesimaHy 
small 5k, a- The correlation length and similar quantities can be computed for arbitrary 
5 k a and X\, and the interplay between non- universal and universal behaviour can be 
studied by approaching the critical temperature (5k\ — > 0). 

In this paper we discuss these ideas for a specific model of two coupled fields with 
discrete symmetry. The parameter x\ will be a measure of the strength of the coupling 
between the two fields. For the present investigation our main purpose is not high accuracy 
for the universal quantities but rather a demonstration that numerical solutions of flow 
equations can give a realistic picture. We, therefore, limit ourselves to a relatively rough 
truncation of r^. More involved truncations would give better quantitative accuracy 
but also obscure the presentation of the basic mechanisms and complicate the numerical 
work. Such improvements will become relevant once a given experimental setting with 
given form of Ta is investigated, with emphasis on quantitative accuracy. They are also 
needed if one attempts to compute universal quantities such as crossover exponents with 
an accuracy comparable with other existing methods. 

2. The evolution equation for the potential 

We consider a theory of two real scalar fields Xa (a = 1,2), in (^-dimensional Euclidean 
space, with an action S[x] invariant under the discrete symmetry (1 «-> —1, 2 <-> —2, 1 «-> 
2). We specify the action together with some ultraviolet momentum cutoff A, so that the 
theory is properly regulated. We add to the kinetic term an infrared regulating piece [[| 

^S=\j -^ d R k {q)xl{q)X a {ql (2.1) 

where x a {q) are the Fourier modes of the scalar fields. The function R^ is employed in 
order to prevent the propagation of modes with characteristic momenta q 2 < k 2 . This 
can be achieved, for example, by the choice 

Rk{q)= ^W)' with /fc2(g)=exp (-^)- (2 - 2) 

The quantity Zk is an appropriate wave function renormalization whose precise defini- 
tion will be given below. The modes with q 2 ^> k 2 are unaffected by the infrared cutoff, 
while the low frequency modes with q 2 <C k 2 are cut off, as Rj~ acts like a mass term: 
lim ? 2^ Rk{q) = Zkk 2 . Through the usual Legendre transformation we obtain the generat- 
ing functional for the 1PI Green functions r^f^ 61 ], where <p a is the expectation value of the 
field x a m the presence of sources. The use of the modified propagator for the calculation 
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of T k results in the effective integration of only the fluctuations with q 2 > k 2 . Finally, the 
effective average action is obtained by subtracting the infrared cutoff piece 

t, r, al * r,a, 1 f A 



/ Ttt^qmqWiq). (2-3) 



J 2 7 (2tt) c 

For equal to the ultraviolet cutoff A, T k becomes equal to the classical action S (no 
effective integration of modes takes place), while for k — > it tends towards the effective 
action F (all the modes are included), which is the generating functional of the 1PI Green 
functions computed from S without infrared cutoff. For intermediate values of k the 
effective average action realizes the concept of a coarse grained free energy ||. 

The interpolation of T k between the classical and the effective action makes it a very 
useful field theoretical tool. The means for practical calculations is provided by an exact 
flow equation, which describes the response of the effective average action to variations 
of the infrared cutoff (t = \n(k/A)) [g 

= |(ri 2) [0] +R k r 1 ^R k ^ . (2.4) 

(2) 

Here 1^ is the second functional derivative of the effective average action with respect 
to <p a . 

For the solution of eq. (|2.4| ) an efficient truncation scheme has to be developed. 
The form of the effective average action is constrained by the (1 <-> —1,2 <-> —2,1 «-> 
2) symmetry. However, there is still an infinite number of invariants to be considered. 
Throughout this paper we work with an approximation which keeps only a standard 
kinetic term in the effective average action 

v k = j d d x jf4( Pl , P2 ) + l -z k d^ a d,r) , (2-5) 

and neglect all invariants which involve more derivatives of the fields. We have used the 



definition p a = ^(p 2 . The parameter appearing in eq. ( |2.2|) can be identified with the wave 
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function renormalization Z k . The justification for our approximation lies in the smallness 
of the anomalous dimension, which is expected to be rj ~ 0.04 for the three-dimensional 
theory. We estimate the corrections arising from the inclusion of higher derivative terms 
to be of the same order as 77. In order to obtain an evolution equation for the effective 



average potential U k from eq. (|2 .4| ) , we have to expand around a constant background field 
configuration (so that the derivative terms in the parametrization ( |2.5|) do not contribute 
to the l.h.s. of eq. Q). Equation (EQ) then gives [EL 



d If d d q ( 1 1 \ d 

di Uk{phP2) = 2 J (2^F {P(q 2 ) + M 2 + P(q 2 ) + M| J dl Rkiq) - (2 - 6) 

P{q 2 ) results from the combination of the standard kinetic contribution Z k q 2 and the 
regulating term R k into an effective inverse propagator (for massless fields) 

P(q 2 ) = Z k q 2 + R k = Zk f (2.7) 

1 — Jh \Q) 
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with f%(q) given in eq. ( |2.2j ); Mf 2 are the eigenvalues of the mass matrix at the point 
(Pi,p 2 )l|: 



M i 2 = -x\U 1 + U 2 + 2U u pi + 2f/ 22 p 2 ± {Ui - U 2 + 2U u pi - 2U 22 p 2 ) 2 + l6U 2 2 p lP2 



(2.8) 

and we have introduced the notation U\ = dUk/dpi, U\ 2 = d 2 Uk/dpidp 2 , etc. 

The fixed point structure of the theory is more easily identified if we use the dimen- 
sionless renormalized parameters 

p a = Z k k 2 ~ d p a 



m 2 



Z^k~ 2 Ml 



(2.9) 



The rescaled mass eigenvalues m\ 2 are related to Uk,pi, p 2 through expressions analogous 
to eqs. (|2.8|), with u\ = duk/dpi, etc. The evolution equation (|2.6|) can now be written 



in the scale-independent form 

d 



dt 



u(pi, h) = ~du +(d-2 + r])(p lUl + ~p 2 u 2 ) - v d L d Q {m{) - v d L%(m 2 2 ), (2.10) 



with 



2 d+ \i r (i 



The functions Lq(w) are given by 



(2.11) 



L d (w) = -2 



dyy z 



exp(-y) 
[1 - exp(-y)]' 



We also employ the functions 



d 



1 8 

n ow 



1 - exp(-y) 



+ w 



for n > 1 



(2.12) 



(2.13) 



in the following. These functions are negative for all values of w. Also |L^(w)| are mono- 
tonically decreasing for increasing w and introduce threshold behaviour in the evolution. 
For large values of m 2 ~ M 2 / k 2 the last two terms in eq. ( |2.10[ ) vanish and the evolution 
of Uk stops. The above functions have been extensively discussed in refs. || |7|, [T^|. As 
they have no simple analytical form, we use numerical fits in the following. 

The evolution of the wave function renormalization Zk, which gives the anomalous 
dimension 

d(lnZ k )/dt = - V , (2.14) 
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can be computed from the exact flow equation ( |2.4[ ) by expanding around a background 
field configuration with a small momentum dependence . We work in three dimensions 
and use the approximate value rj = 0.04. Within the accuracy of our results, this value is 
consistent with the values of r\ in both the Ising and Heisenberg models, which are relevant 
near the fixed points. We take into account threshold effects in the evolution of and 
the vanishing of rj in the symmetric regime (when the minimum of is at the origin) 
by using rj = 0.04 as long as the theory is in the regime with spontaneous symmetry 
breaking and the mass eigenvalues are in the interval 0.01 < fh 2 a < 100. Otherwise we set 
f] = 0. We also point out that we have neglected a term proportional to 77 in the r.h.s. of 
( [2.10 ), which comes from the t-derivative acting on inside Rk in the r.h.s. of eq. (|2.6|) . 
This is consistent with our approximate treatment of the wave function renormalization 
effects and is justified by the smallness of rj. 



3. The fixed point structure 



From this point on we concentrate on the three-dimensional theory. Equation ( |2.10 ), with 



d — 3, is the master equation for our investigation. It is a non-linear partial differential 
equation for three independent variables. A numerical solution is possible (see refs. |I| [□]] 



for the solution of the evolution equation for the 0(iV)-symmetric scalar theory). However, 
the presence of three independent variables complicates the numerical treatment. For 
this reason we resort to a different approximation scheme. We first introduce a Taylor 
expansion of Uk{pi, P2) around a certain point (usually the minimum of the potential or 
the origin). This turns eq. ( |2.6|) into an infinite system of ordinary (coupled) differential 
equations for the /c-dependence of the point around which the expansion is performed and 
the derivatives of the potential, with independent variable t = ln(A;/A). We solve this 
system approximately by truncating at a finite number of derivatives. This approach has 
been used in the past for the study of the 0(iV)-symmetric scalar theory. It has provided 
a full, detailed picture of the high-temperature phase transition for the four- dimensional 
theory [1^, [7|, |13|, [14j , with accurate determination (at the few % level) of such non-trivial 



quantities as the critical exponents 0j. An estimate |16| of the residual errors for high 
level truncations indicates that they are smaller than the uncertainties introduced by the 
imprecise treatment of the wave function renormalization effects. This has been confirmed 
by comparison with the numerical solution of the untruncated evolution equation for the 
potential |T(J [□]]. The same approach has been used in ref. [[J for the study of the 



high-temperature phase transition in a four- dimensional theory of two scalars with a 
(1 «-> —1,2 <-> —2, 1 2) symmetry. The truncation employed there kept up to second 
derivatives of the potential Un, Uu, Uzi- The present work is a detailed investigation 
of the crossover phenomena in the effective three-dimensional theory near the critical 
temperature. We use an extended truncation which keeps up to the third derivatives of 
the potential Um, Uvn-, etc. We also note that this theory belongs in the same universality 
class as the statistical two-component spin systems with cubic anisotropy, whose universal 
behaviour has been studied in p|-p|. 
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Let us first parametrize the potential by its derivatives at the origin (S regime) 

m 2 {k) = E7i(0) = U 2 (0) 

\(k) = U n (0) = U 22 (0) g(k) = U 12 (0) 

Mk) = tfm(O) = U 222 (0) v B {k) = £7 112 (0) = y 122 (0), (3.1) 

where the equality of derivatives is imposed by the discete symmetry of the theory. Our 
approximation amounts to using, for the potential, the ansatz 

™2/„ , „ \ , ^ / 2 i „2\ | -„ „ | / 3 . 3 



£4(pi, pa) = ™ 2 (pi + p 2 ) + - (p 2 + p 2 .) + w 2 + y (PI + P2 j + y (P1P2 + P1P2 j • (3-2) 

Evolution equations for the above parameters are obtained by taking derivatives of eq. (|2.6|) 
with respect to pi >2 . It is convenient to define the dimensionless couplings 

m 2 (k) =Z k l kr 2 m 2 (k) = AT 2 m 2 (k) 

A(jfc) =Zf k' 1 X(k) = k~ l \{k) g{k) = Zf AT 1 g(k) = AT 1 g(k) 
u A {k) =Zf u A (k) v B (k) = Zf u B (k) (3.3) 



The dimensionless potential Uk(pi,p 2 ) defined in eqs. (|2.9|) can be expressed in complete 
analogy to eq. ( |3.2|) in terms of these dimensionless quantities. The evolution equations 
for the above quantites have a scale-invariant form, in the sense that the r.h.s. does not 
explicitly involve a dependence on k 



--{-2 + r])m 2 + v 3 (3\ + g)Ll(m 2 ) (3.4) 

= (_1 + 2r] )\ - v 3 (9\ 2 + f)L\{m 2 ) + v 3 (5u A + v B )L\{m 2 ) (3.5) 

= (_1 + 2r] )g - v 3 {6\g + Ag 2 )L 3 2 (m 2 ) + Qv 3 v B L\{m 2 ) (3.6) 

~-3r]u A + 2v 3 (27\ 3 + g 3 )L 3 (m 2 ) - 3v 3 (lb\v A + gu B )L 3 (m 2 ) (3.7) 



dm 2 
~dT 

d~X 

~dt 

dg_ 

dt 
dv A 
dt 

=3r]u B + 2v 3 (9X 2 g + 15Xg 2 + A~g 3 )L 3 (m 2 ) - v 3 (2l\u B + h~gu A + 22gis B )L 3 (m 2 ), 

(3.8) 

where v 3 = 1/87T 2 . They automatically preserve the discrete symmetry. 

We are interested in the fixed points of the last set of equations. First there is the 
ultraviolet attractive Gaussian fixed point (G) with m 2 , A, g, u A , v B = 0. There are also 
three fixed points with at least one infrared attractive direction ||17|1-||19|1. They all appear 
for m 2 < 0, A, g, u A , v B > 0, which means that the minimum of Uk(p~i, p 2 ) is away from the 
origin. (The exact values are not important since the discussion in this section concen- 
trates on the qualitative behaviour only.) For their identification we use their standard 
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names in statistical physics [P^]-JI^[: 

a) The Heisenberg fixed point (H) has g = A, z/g = z/4 and corresponds to a theory with 
symmetry increased to 0(2). The potential Uk(pi,p2) has a series of degenerate minima 
at pio + P20 = Po- 

b) The Ising fixed point (I) has g = Vb — and corresponds to two disconnected Z 2 - 
symmetric theories. The potential Uk(p\,p2) has minima at p 10 = p 2 o = Po/2 (for positive 
and negative values of 4>\,2)- 

c) The Cubic fixed point (C) has g = 3A, v B = 5ua and corresponds to two disconnected 
^-symmetric theories, if the fields are redefined according to 

# = -^(0i + 2 ), ^=-^(0x-0 2 ). (3.9) 

The potential Uk(pi, P2) has minima at pio = po, P20 = and P20 = Po, Pio = 0. 
All these points are infrared unstable in the m 2 direction and are located on a critical 
surface m 2 r = m? cr {\,g,VA- l VB) < 0. Solutions of the evolution equations, which start 
above the critical surface, with m 2 > m 2 r , flow towards the region of positive m 2 for 
t — > — 00, and correspond to theories in the symmetric phase. Solutions with ifi 2 < m 2 r 
flow deep into the region of negative m 2 and correspond to theories in the phase with 
spontaneous symmetry breaking. 

The relative stability of the fixed points on the critical surface determines which one 
governs the dynamics of the phase transition in the immediate vicinity of the critical 
temperature (or mass). It has been discussed in ref. 0, in a truncation that neglects 
derivatives of the potential higher than the second (and corresponds to va-,^b = 0). 
All three fixed points are attractive in the A direction. However, the Ising and Cubic 
fixed points are repulsive in the g direction, while the Heisenberg fixed point is totally 
attractive. There are four disconnected regions: 

a) 3A > g > A : The trajectories flow away from the Cubic towards the Heisenberg fixed 
point. 

b) A > g > : The trajectories flow away from the Ising towards the Heisenberg fixed 
point. 

For parameters in the above two regions we expect second-order phase transitions, with 
critical dynamics governed by the three fixed points. 

c) g > 3A : The trajectories flow away from the Cubic fixed point into a region of large g 
and small A. Eventually A turns negative at a finite value of k. At this point an instability 
arises which signals the presence of a first-order phase transition. 

d) g < : The trajectories flow away from the Ising fixed point and cross the line g = —A 
at a finite k. This again implies the presence of an instability whose true nature is related 
to a first-order phase transition. 

Flows that start on the lines g = 0, A, 3A remain on these lines. No trajectories exist 
which connect the four regions g > 3\, 3\ > g > A, A > g > 0, g < 0. 

We are interested in studying in detail the crossover from the Ising to the Heisenberg 
fixed point in the region A > g > 0. The region 3 A > g > A can be mapped on the region 
A > g > through a redefinition of the fields according to eqs. (|3.9| ). This shows that 
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the Ising and Cubic fixed points lead to an identical universal behaviour (and therefore 
to identical universal quantities, such as critical exponents and crossover curves). 



4. The crossover 



We would like to study the crossover from the Ising to the Heisenberg fixed point. In the 
vicinity of these fixed points the potential M fc (p 1 ,p 2 ) has minima at p 10 = p 2 o = p~o/2 > 
(for positive and negative values of 0i, 2 )- In order to improve our quantitative accuracy 
in this regime (which we call the M regime), it is convenient to parametrize the potential 
in terms of its derivatives at the minimum. For this reason we define the parameters 



Hk) = Unfa) 
V A (k) = Uni(po) 



U 2 2(po) 
= U 2 22(po) 



g(k) = Mpo) 
u B {k) = U U2 (po) 



U 12 2 (pO 



(4.1) 



where pio = P20 = Po/2 is the minimum of Uk(pi, P2), related to the minimum of Wfc(pi, p 2 ) 
through the first of eqs. (|2.9| ). The equality of derivatives is imposed by the (1 <-> —1,2 <-> 
—2, 1 <-* 2) symmetry of the theory. The first derivatives U\{po) = U 2 {po) are zero f° r 
non-zero po- Our approximation amounts to using for the potential the ansatz 



U k {pi,p 2 ) 



A 
2 

+ 



Pi ~ 



Po 



+ IP2- 



Pi 



Po 
2 



Po 



+ IP2- 



+ g(pi- 



Po 



P2 



Po 
2 



Po 



+ — — ■ 

2 



Po 



Po 



Po 



Po 

P2 "T 



(4.21 



We define dimensionless couplings according to eqs. (|3.3|) for the derivatives and 

K (k) = Zkk^p^k) (4.3) 

for the minimum. The evolution equations for these couplings are lengthy and can be 
found in the appendix. They automatically preserve the discrete symmetry. The various 
fixed points are now located on a critical surface n 2 cr = ^(A, g, Ua, Vb), which is unstable 
in the k direction and separates the phase with spontaneous symmetry breaking from the 
symmetric one. Solutions of the evolution equation starting with k > n cr flow towards 
the region of large k and correspond to the phase with symmetry breaking. Eventually 
the minimum of the potential settles down at a finite value 



Poi? = limiTfcPofA;) = lim &/«(&). 

fe^O fc^O 



(4.4) 



Solutions starting with k < n cr lead to n(k) = at some finite k s . From this point 
on, the system is in the symmetric regime (S regime) and we use a parametrization of 
the potential according to eqs. ( |3.1|) and (|3.2|) , the rescaled variables defined in eqs. 
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Q3.3| ), and the evolution equations ( |3.4j )-( pUf ). Finally, if the evolution starts on the 



critical surface the flows never lead out of it. Instead they eventually approach the most 
attractive fixed point on the critical surface. This means that for k — > the renormalized 
minimum of the potential is given by eq. (|4.4j) with finite k. Therefore, poi? = and the 
critical surface corresponds to a second-order phase transition. 

We numerically integrate the evolution equations by starting at some high momentum 
scale (ultraviolet cutoff) k = A (t — 0) where the potential U\ is equal to the classical 
potential V. We consider U\ quartic in the fields, and therefore given by eq. ( f4.2| ) with 
Pa (A) = Pb(A) = 0. For the wave function renormalization we start with Z\ — 1. We 
fine-tune p (A)/A = k(A), so that \5k\\ = \k(A) — K cr \ < 1, and the system is very close 
to the critical surface. The subsequent evolution moves the parameters along the critical 
surface and close to various fixed points depending on their relative stability. We define 
the quantity 

in order to parametrize the distance from the Heisenberg fixed point. We have xi — — 1 at 
the Ising fixed point and xh = at the Heisenberg one. We can arrange for the evolution 
to approach first the Ising fixed point by selecting 6x\ = x(A) + 1 <C 1. Then x(k) 
remains very close to its initial value while the running couplings approach their values 
at the Ising fixed point 

Ki = 8.57 x 1(T 2 , A/ = 8.05, gj = 0, (v A )i = 75.8, {y B )i = 0. (4.6) 

Near the fixed point the evolution slows down, as the /^-functions are almost zero. How- 
ever, the Ising fixed point has two unstable directions and the system will eventually 
move away from it. One of these directions leads away from the critical surface, while the 
other is directed towards the Heisenberg fixed point. By selecting the starting point of 
the evolution so that 5x^ 3> we can arrange for the flow to stay close to the critical 

surface while moving from one fixed point to the other. Eventually the running couplings 
approach their values at the Heisenberg fixed point 

k h = 6.02 x 10~ 2 , A H = 7.32, g H = 7.32, (u A ) H = 56.4, (v B ) H = 56.4. (4.7) 

The trajectory in coupling space connecting the two fixed points is depicted in fig. 1, 
while the /3-functions for the various couplings along the trajectory are plotted in fig. 2. 
It is clear that the evolution is very slow near the fixed points where the /3-functions are 
zero and fast in the middle of the trajectory. As long as \Sk,\\/Sxa «C 1 the system stays 
close to the critical surface for a sufficiently long "time" t for the Heisenberg fixed point 
to be reached. However, if the ratio \5k\\/5x\ is increased the system leaves the critical 
surface at some earlier point of the evolution along the critical trajectory of fig. 1, and 
moves towards the phase with symmetry breaking or towards the symmetric one. For 
5k^ < it runs into the S regime at some finite k s and eventually settles down in the 
symmetric phase with renormalized parameters [0, [J 

m 2 R =limm 2 (/c) = lim k 2 m 2 {k) 

k *0 k~ ^0 
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9r 

Xr — lim X(k) = lim kX(k) gR = lim g{k) = lim kg{k) xr = 1 

k— >0 fc^O k^O k^O Xr 

v AR = lim u A (k) v B r = lim u B (k). (4.8) 

fc^O k^O 

For 8k\ > the system settles down in the phase with symmetry breaking, with the 
minimum of the potential given by eq. ( |4.4|) and the various couplings at the minimum 
defined in analogy with eqs. (f4.8| ) []. 

The singular behaviour of the theory in the critical region near the phase transition 
can be parametrized in terms of critical exponents. For example, the renormalized mass 
in the symmetric phase can be expressed as 

m\ oc \5k a \ 2v . (4.9) 

The parameter v can be computed as the slope of the function lnm^ = / (In |5/ca|)- If 
the system is in the vicinity of a fixed point, v is constant and can be identified with the 
corresponding critical exponent. Moreover, this exponent can be related [[T], ^|, |2D| to the 
most negative eigenvalue Ai of the matrix /3y = dfii/dcij (with a,j = k, X, g, ua, vb), which 
describes the growth of small perturbations 5n = k — kj p , 5X = X — Xf p , etc., around the 
fixed point, i.e. 

v = ~. (4.10) 

In fig. 3 we plot the three smallest eigenvalues of the matrix flij for small perturbations 
5k = k — Kt r , SX = X — X tr , etc., around the trajectory of fig. 1. For this calculation we 
employ the /^-functions given in the appendix. We observe that the smallest eigenvalue 
remains negative along the whole trajectory. It corresponds to the unstable direction 
away from the critical surface. From its values at x = —1 and x — 0, making use of 
eq. (|4.10|) , we can infer the critical exponent v of the Ising and Heisenberg models 

vi =0.700 (4.11) 
Ujj =0.733. (4.12) 

The above values deviate at the 10% level from the more accurate results obtained through 
the e-expansion or summed perturbation series (i/j = 0.631(2), vn = 0.670(3) PH), or from 
results obtained with the same method and more extended truncations M. This is not 
surprising, as our very rough treatment of the wave function renormalization generates 
uncertainties of the order of the anomalous dimension t], as we mentioned at the end of 
section 2. We observe that the difference Ujj — uj = 0.033 compares rather well with the 
high precision estimate 0.039. We have verified the above values for the critical exponents 
by the explicit computation of the function lnm^ = / (In |5ka|). For this we integrated 
the evolution equations for theories with \8k\\,8xa <C 1, as we discussed in the previous 
paragraphs. For \5k\\ 3> 8x\ the evolution approaches the Ising fixed point and slows 

5 In the vicinity of the Heisenberg fixed point the presence of a massless Goldstone mode results in 
zero values for A^,, gu in the limit k — » 0. Quartic couplings can be defined at non-zero values of k [12l m|. 
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down for a long "time" t. As a result the system loses memory of the initial conditions 
of the running (the classical parameters of the theory). Eventually the system leaves the 
critical surface in the vicinity of the Ising fixed point. The function lnm^ = / (In |<5/ca|) 
has a constant slope, which reproduces the value of eq. (|4.11[) for the exponent v. For 
theories with a smaller |5ka|/&£a ratio the evolution first approaches the Ising fixed point 
and then deviates from it along the critical surface. Finally it leaves the critical surface 
at some point along the trajectory of fig. 1 connecting the two fixed points. The function 
lnm^j = / (In |5«a|) now has a variable slope, since the growth of deviations from the 
critical surface depends on the average of the eigenvalue Ai along the part of the trajectory 
covered by the evolution. Only for |5/«a| •C &ca does the slope become constant again. For 
these theories the flow first moves from the Ising to the Heisenberg fixed point, and then 
spends a long "time" t in its vicinity. During this "time" it loses memory of the previous 
evolution, before eventually leaving the critical surface. The calculation of the exponent 
v reproduces the value of eq. (|4.12|) . The error in the above values of v originating in 
uncertainties of the numerical integration is less that 1% and, therefore, much smaller 
than the 10% error resulting from the truncations of the evolution equation. 

We now turn to the second smallest eigenvalue of fig. 3, which changes sign along 
the trajectory connecting the two fixed points. Its negative sign at x = —1 indicates 
that the Ising fixed point has two unstable directions, while its positive value at x = 
indicates that the Heisenberg fixed point has only one. The eigenvalue A2 corresponds 
to the direction along the critical surface, which leads from the Ising to the more stable 
Heisenberg fixed point. Its value at x = —1 determines the growth of perturbations 
that couple the two disconnected ^-symmetric theories of the Ising fixed point. For the 
renormalized theory these perturbations are parametrized in terms of the dimensionless 
quantity xr + 1 = Qr/Xr- A crossover exponent ip can be defined through the relation 

x R + 1 oc \5k a \- v . (4.13) 

In the vicinity of the Ising fixed point this exponent is related to the two smallest eigen- 
values according to |l|, p|, |20| 



A? 

<p = = 0.032. (4.14) 
Ai 

Alternatively ip can be computed from the slope of the function \h{xr + 1) = / (In |5«a|)- 
The two results are in agreement. Comparison with more accurate results obtained 
through the e-expansion {ip = 0.110(5), |T7], [0|) indicates that our approximations result 
in limited accuracy for such a small quantity. In fact, the eigenvalue | A2 1 is accidentally 
very small on the scale of the other eigenvalues |Ai| and | A3 1 (cf. fig. 3). A modest trun- 
cation error can easily change it by a factor 3 or 4. It is not very difficult to increase the 
truncation level in order to reach a few per cent accuracy for v J7| and correspondingly 
for ip. (We guess that a truncation at the level of ref. J7] should yield eigenvalues Aj with 
an accuracy of ~ 0.05 and, therefore, an uncertainty in the ip value ~ 0.03.) A precise 
estimate of <p is, however, not the purpose of this paper. As mentioned in the introduction 
we rather want to demonstrate that our method is capable of covering the whole crossover 
region between the two fixed points. 
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In fig. 4 we plot the difference of the "effective" exponent v = dlnrriR/dln \ 8ka\ from 
its value at the Ising fixed point as a function of x(A). For this we integrate numerically 
the evolution equations by starting at the ultraviolet cutoff k — A (i = 0) where the 
potential Ua is equal to the classical potential V. We consider U\ quartic in the fields, 
with A (A) = 5 and fine-tune k(A) so that the flows approach the critical surface. By 
varying \8k\\ we obtain theories with various values of the renormalized mass. In fig. 



4 we also plot the renormalized parameter xr defined in eqs. (|4.8|) . The curves (a) 
correspond to theories with renormalized mass ttir/A ~ 10~ 4 - 10~ 3 , while the curves (b) 
correspond to theories with renormalized mass trr/ A ~ 10~ 5 - 10~ 4 . The presence of 
small kinks in these curves gives a feeling of the typical errors resulting from the numerical 
integration - a less than 1% effect. It is clear that the evolution results in values for x 
in the limit k — > which are closer to the Heisenberg fixed point than the "initial value" 
x(A). Both curves for xr are above the dotted line, which corresponds to xr = x(A). 
Theories with smaller renormalized mass stay close to the critical surface for longer "time" 
t, so that the attraction of the Heisenberg fixed point is more pronounced. For given x(A) 
this results in xr being closer to zero for curve (b) than for curve (a). Similarly the value 
of v for curve (b) is closer to the Heisenberg fixed point value than for curve (a). We 
have also included in fig. 4 the line corresponding to —\{ 1 (x) (dashed-dotted line), with 
Ai the most negative eigenvalue of the stability matrix plotted in fig. 3. For sufficiently 
small ttir/A the flows quickly approach the critical trajectory of fig. 1 at some point 
x ~ x(A). Then the running follows the critical trajectory towards the Heisenberg fixed 
point, until eventually it moves away from the critical surface towards the symmetric 
regime. The resulting exponent v corresponds to an average over the values of —\{ x [x) 
for the part of the flow close to the critical trajectory, along with some non-universal 
contributions from the initial and final running. For sufficiently small m^/A (such as for 
curve (b)) the latter contributions become insignificant. The resulting critical exponent 
v is larger than the value of — A^f 1 at the point x(A), because of the running towards 
the Heisenberg fixed point. For this reason the curve (b) for v lies above the dash- 
dotted line. For the curve (a) the non-universal initial and final running give significant 
contributions, which contaminate the universal part along the critical trajectory. The 
resulting "effective" exponent v has a complicated dependence on x(A). It should be 
pointed out that in most experimental situations the observed values of tur/ 'A are not 
sufficiently small for the non- universal physics to be insignificant. It seems difficult to 
perform investigations of crossover based only on the value of the critical exponents. This 
would require to disentangle the crossover from the non-universal behaviour. Figure 4 
is an explicit demonstration that our method can effectively deal with such situations. 
Starting from a given set of values for the parameters of the theory at the short distance 
scale A, both the universal and non- universal long distance physics can be determined 
through the integration of the evolution equations. 



12 



5. Conclusions 



We have demonstrated that the flow of the effective average action T k can be used for 
a description of crossover phenomena over the whole crossover region. We have com- 
puted both universal quantities - the eigenvalues \(x) of the stability matrix of small 
fluctuations along the crossover curve parametrized by x - as well as actual "effective" 
exponents (see fig. 4) for given values of m R /A. Even though an exact non-perturbative 
flow equation for T k is available, practical calculations are feasible only for truncations of 
rV These approximations have not been pushed to a high level of accuracy in the present 
paper. High precision calculations within the e-expansion have determined the crossover 
exponent tp with very satisfactory accuracy in the past. However, for many practical 
applications the actually observed critical behaviour in crossover situations is not only 
determined by the universal critical behaviour. It is often a mixture of the flow towards 
the fixed points and the subsequent crossover flow from one fixed point to another. In 
these situations of "imperfect scaling", many standard methods, as for example the e- 
expansion, cannot be applied easily. The flow of the effective average action is directly 
formulated in three dimensions and includes the relevant degrees of freedom at every scale 
k (not only for k in the scaling region). It can be used both in the symmetric phase and 
in the phase with spontaneous symmetry breaking, where the existence of massless Gold- 
stone bosons hinders the application of many standard methods. For a given initial form 
of the "microscopic action" T A , we can solve these equations numerically and compute 
the quantities of interest. We hope that the present method will develop in the future 
into a useful practical tool for quantitative investigations of crossover problems. 
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Appendix: The evolution equations in the M regime 

It is easier to derive the evolution equations for the dimensionless derivatives in eqs. ( |3.3|) 
and the rescaled minimum of the potential of eqs. ( f4.3|) by making use of the evolution 
equation ( 2.10|) for the rescaled potential Uk{p\,p~2) in three dimensions 



d_ 

Of 



u(pi, fa) = - 3n + (1 + r})(jp\u\ + p 2 u 2 ) + C 



( = -v 3 Ll{ml)-v 3 L 3 {ml), 
where the rescaled mass eigenvalues are given by 



(A.l) 



m i iifti, fa) =o ( M i + u 2 + 2u n pi + 2u 22 p2 ± -4) 



.4 



[ux -u 2 + 2ttnpi - 2u 22 p 2 ) + 16u 12 p 1 p 2 



(A.2) 



The evolution of n(k) is obtained by taking the total t-derivative of the minimization 
conditions duu/dpi\ K = duk/dp 2 \ K = 0. We obtain 



eft" {+V) X + gdfa 



For the derivatives we find 

dX 
~di ' 



-l + 2r])X + ^ 



dg 
dt 



^=(-1+2^ 



dpi 

d 2 ( 



v A + v B d( 



dpidp 2 



X + g dpi 
2v B d( 



X + g dpi 



dUA 

dt 

dv B 
dt 



-3r]V B + 



dpi 



d~p\d~p 2 



The necessary p-derivatives of ( are given by 
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d 2 C 
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v 3 



^3 
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dm 2 , 
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dp 



X2 
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dpi dp 2 



dpi dp 2 
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(A.4) 



14 



+ ^3 



dp 1 dp 2 



dp 1 dp 2 



dpi 



2v, 



dmf 
dpi 



dm 2 , 
dpi 



Ll{rhl) 



- 3% 



+ ^3 



dm\ d 2 m\ , , , _ , dm 2 d 2 m 2 
dpi dpi 



2 r3/ ~ 2\ 



dpi dp( 



d 3 m 2 



d 3 



~ 2 

m 



dpi 



dp 



d 3 ( 
dpidp2 



2v, 



dmi 
dpi 



d~p 2 Ls(mi) + 



dm? 
dpi 



dp 2 



-2v 3 



dmi d 2 mi T?J( ~ 2 \ dm\ d 2 m^ 



^3 



+ ^3 



2 L\{m\) 



dpi dpidp 1 
h\ dm 

dpi dp 2 ~' zv '" u 1 dpi dp 2 

9^ Ll(mi) + ^f^ Ll(m2) | 



9pi dpidp 2 

d 2 mi dmi 3 , „ 2*. d 2 m 2 dm 2 
~i^2 — ~x^L 2 (m l 



Ll(ml) 



(A.5) 



At the minimum the mass eigenvalues of eqs. (|A.2|) and their derivatives are given by 
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Figures 



Fig. 1 : The trajectory on the critical surface which leads from the Ising to the Heisenberg 
fixed point. 

Fig. 2 : The /3-functions for the various couplings along the trajectory of fig. 1. 

Fig. 3 : The three smallest eigenvalues for perturbations around the trajectory of fig. I. 

Fig. 4 : The difference f3 — (3i and the renormalized parameter xr as a function of the 
value of x at the cutoff. Lines (a) correspond to theories with renormalized mass 
thr/A ~ 10~ 4 - 10~ 3 , while lines (b) correspond to theories with renormalized mass 
m R /A ~ 10~ 5 - KT 4 . 
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